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Abstract. The formation of the optical polaron and bipolaron in two-dimensional 
(2D) systems are studied in the intermediate electron-phonon coupling regime. The 
total energies of 2D polaron and bipolaron are calculated by using the Buimistrov- 
Pekar method of canonical transformations and analyzed in the weak, intermediate and 
strong coupling regimes. It is shown that the electron-phonon correlation significantly 
reduces the total energy of 2D polaron in comparison with the energy of the strong- 
coupling (adiabatic) polaron. A charge carrier in polar crystals remains localized in a 
2D potential well when the electron-phonon coupling constant a is greater than the 
critical value a c ~ 2.94, which is much lower than a critical value of the electron- 
phonon coupling constant a for a 3D system. The critical values of the electron- 
phonon coupling constant a and the parameter of Coulomb repulsion between two 
carriers (3 = 1/(1 — £oo/£o) (where and £0 are the high frequency and static 
dielectric constants, respectively), which determine the bipolaron stability region, are 
numerically calculated. The obtained results are compared with the ones obtained 
by using the Feynman path integral method and the modified Lee-Low-Pines unitary 
transformation method. 
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1. Introduction 

As is well known, charge carriers (electrons and holes) is ionic crystals interact with 
polar optical phonons and they are self-trapped at their sufficiently strong coupling to 
the phonons with the formation of optical polarons (called as the Pekar or Frohlich 
polarons) [H EJ [3] . One distinguishes three distinct regimes of electron-phonon coupling 
[2]: (i) the weak-coupling regime describes the correlated motions of the lattice atoms 
and the quasi-free charge carriers which remain in their initial extended state, (ii) the 
intermediate-coupling regime characterizes the self-trapping of a charge carrier which 
is bound within a potential well produced by the polarization of the lattice in the 
presence of the carrier and follows the atomic motions, and (iii) the strong-coupling 
regime describes the other condition of self-trapping under which the lattice atoms 
no longer follow the charge carrier motion and the self-trapping of carriers is usually 
treated within the adiabatic approximation (i.e. lattice atoms remain at their fixed 
positions). Under certain conditions, two charge carriers interacting with the lattice 
vibrations and with each other can form a bound state of two carriers in polar materials 
within a common self-trapping well. Since the attractive interaction of the electron- 
phonon coupling in these systems is strong enough to overcome the Coulomb repulsion 
between two carriers. The self-trapped state of the pair of charge carriers is termed 
a bipolaron. In the last few decades, the bipolaron problem has been a focus of 
attention due to its importance in semiconductor technology and in the bipolaronic 
mechanism of superconductivity (see Refs. [H 13, EJ d El [9], UTIl HH H2])- After the 
discovery of the layered high-T c cuprate superconductors, the study of bipolarons has 
attracted the revived interest because some researchers believe that the bipolaron is 
one of the possible candidate for explaining high-T c superconductivity. The mechanism 
of high-T c superconductivity of large bipolarons was proposed by Emin and Hillary 
[7j. This mechanism is based on the Bose-Einstein condensation of bipolarons as 
discussed by Schafroth [13]. Another mechanism of high-T c superconductivity driven 
by the superfluid single particle and pair condensation of large bipolarons and polaron 
Cooper pairs was proposed in [TJ1E5]. The possibility of such a novel superconductivity 
depends on the existence of polaron and bipolarons in the superconducting materials. 
The above electron-phonon coupling regimes are characterized by the dimensionless 
Frohlich coupling constant a. The possible ranges of the values of a characterizing 
the weak, intermediate and strong coupling regimes depend on the type of self-trapped 
charge carriers and the dimensionality of the system. The polaron and bipolaron ground 
state in the strong coupling limit (i.e. in the adiabatic approximation) have been studied 
by many authors for 3D and 2D systems 0, [7J El EJ EH1 EH ES Ej- However, with 
decreasing a, it is necessary to take into account the electron-phonon correlation which 
reduces the energies of the polaron and bipolaron in comparison with those of the 
adiabatic polaron and bipolaron. Therefore, a quantitative treatment of the (bi)polaron 
problem in the intermediate-coupling regime was also a subject of numerous studies in 
the past three decades (see Refs [61 [8j Ell ESI HI EH]). The ground state energies 
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of the polaron and bipolaron have been calculated by using several approximations, 
such as the Feynman path integral method [HI [9j [101 [HI [20], the Lee- Low-Pines (LLP) 
unitary transformation method [321 [18] and the operator formalism [HI [IE]. According to 
different intermediate coupling treatments, the formation of the polaron and bipolaron 
becomes possible only if the values of a are greater than the certain critical values a c . 
The previous calculations have shown that the (bi)polaron is created more easily in 
2D systems than in 3D ones and the stability region for bipolaron formation is much 
broader in 2D compared with 3D case. So far, calculations of the ground-state 

energies of the (bi)polaron, the values of a c and the stability region of the bipolaron 
are not conclusive. Some of the estimated values of a c differ greatly in magnitude. 
Further, the calculation methods used to study the (bi)polaron problem suffer from the 
certain drawbacks (see Refs. 0,118]). In particular, the Feynman path integral method 
is more accurate for the calculation of the polaron energy pi [161 HB], but it fails to 
describe correctly the bipolaron formation. Moreover, the functional for the bipolaron 
energy has rather a cumbersome form which requires the tedious numerical calculations. 
Therefore, it is expedient to use the another method for the calculation of the polaron 
and bipolaron energies in the intermediate-coupling regime. One of such methods is 
the Buimistrov-Pekar method |21j. In this approach, the electron-phonon correlation is 
taken into account through the displacement amplitude F$(r) of the form: 

F <i(r) = k + 9qexp (-iqr) , ( 1 ) 

where f$ and are the variational parameters, which are determined from the condition 
for minima of the total energy, f is the electron coordinate, if is the phonon wave vector. 

The second term in Eq. ([Q) takes into account the electron-phonon correlation 
effect that leads to the reduction of the energy of a 3D adiabatic polaron when a > a c . 
At a < a c the present method leads to the delocalized state of a polaron, and the 
total energy of the polaron is given by E p = afru (where u Q is the longitudinal optical 
(LO) phonon frequency). In the Feynman path integral variational approach [20], the 
derealization of the polaron at a finite value of a does not occur and the lowest polaron 
energy is obtained. However, at a > a c the Buimistrov-Pekar method taking into 
account an important part of the electron-phonon correlation gives also reasonable 
results for the energy of the polaron. The main advantage with this method is that it is 
simple and does not require tedious numerical calculations. Recently, the Buimistrov- 
Pekar method was applied to study both the free bipolaron and the bound bipolaron in 
3D systems [22l l23l [24] . The existence of polarons and bipolarons in high-T c cuprates 
and other materials has been indicated by several experiments [251 [261 [271 1281 129] . The 
Cii0 2 -based layered high-T c materials are believed to be quasi-two-dimensional systems 
and in an intermediate electron-phonon coupling regime. So far, the problems of the 
2D polaron and bipolaron are not studied sufficiently and the possibility of formation 
of intermediate-coupling polarons and bipolarons in 2D systems within the Buimistrov- 
Pekar formalism is not explored. The aim of the present paper is to study the formation 
of the intermediate coupling optical polaron and bipolaron in 2D systems by using 
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the Buimistrov-Pekar method and to expose the important features of this method. 
We calculate the ground-state energies of the intermediate coupling 2D (bi)polarons, 
the critical values of a, the parameter of the Coulomb repulsion between two carriers 
(3 = 1/(1 — Eoo/eq) (where £ ro and e are the high frequency and static dielectric 
constants) for the formation of 2D polaron and bipolaron. We discuss the obtained 
results and compare them with previous ones obtained by using the Feynman path 
integral method and the modified LLP unitary transformation method [16, 30J. 



2. Formation of intermediate-coupling polarons 



The Hamiltonian and variational wave function describing the interacting system of 
electron (or hole) and LO phonons can be written as [221 El 

h 2 
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m* is the effective mass of the carrier before polaron formation, Vg- is the 2D form factor 
of the electron-phonon interaction, L 2 is the size (or surface) of the 2D system, b~t{b q ) 
is the creation (annihilation) operator of the LO phonon with the wave vector q and 
the frequency u , <p{r) is the one-electron wave function represents the phonon part 
of the wave function, |0) is the unperturbed zero-phonon state satisfying b^O) = and 
(0|0) = 1. Averaging the Hamiltonian (JSJ) over $ p ^, we have 
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After substituting Eq.Q into Eq.(j6]), averaging the Himiltonian (jBJ) over (p(r), and 
minimizing the energy E p = (tp\H p \tp) with respect to fg- and g$, we obtain the following 
functional for the total energy of a 2D polaron 
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The latter sum in Eq. j7j) is the correction term arising from the electron-phonon 
correlation caused by the second term in Eq.flT]). The functional (J7J) without this term 
(g$ = 0) determines the energy of the strong-coupling adiabatic polaron. For simplifying 
the calculations, we choose the electron wave function in the form 

<p(r) = Nexp(-5 2 r 2 ), nN 2 = 25 2 J W 2 = exp(-q 2 /A5 2 ) (8) 

Then we obtain the following expression for the energy (in units of hw ) of the 
intermediate coupling 2D polaron: 



o 



r- I [1 - exp(-t 2 )] 2 

E p = 2a* - - 2»a j dt 1 _ eM _ t2)+ ^ H2 , (9) 



where \i = l 8. 

In the weak and strong coupling limits, we obtain from Eq. ([9l) 
~ol for a < a c = 2.94 

E p={ ,2 r ( 10 ) 
1 -far for a — >■ oo 

The Feynman path integral method and the LLP unitary transformation method 
applied to the polaron problem, give the following expressions for the energies of the 
polarons [16]: 
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where if, w and A, 7 are the respective variational parameters. 

The energies of the 2D polarons in the weak and strong coupling limits follow 
directly from Eqs. (11) and (12): 
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Figure 1. The energies of the 2D polarons (in units of tkv ) as a function of a for the 
strong coupling (dotted curve) and weak coupling (dashed line) limits, calculated by 
using Eq. (|10p , for the intermediate coupling limit (solid curve), calculated by using 
Eq.©. Thin solid curve calculated by the Feynman path integral method, Eq. (fTTj) . 



These results can be compared with Eqs. ([9]) and (fit)]) . As can be seen from Eqs. (JTDjl . 
(JT3l) and (JT4l) . all three methods in the strong-coupling limit give the same result. For 
compairing, in Fig.l, the energies of 2D polarons (in units of hu) ) obtained from the 
Eqs. ([9]) and ffTT]) are plotted as a function of a. As seen from Fig. 1, the energy of 
the 2D polaron calculated by taking into account the electron-phonon correlation within 
the Buimistrov-Pekar approximation is lower than that of the 2D adiabatic polaron. At 
a > a c ~ 2.94, a bound state of a carrier exists in 2D systems and the polaron remains 
localized in the 2D polarization well. 



3. Formation of the intermediate-coupling bipolarons 



In the bipolaron problem the total energy of the interacting system of two carriers 
and LO phonons is determined. The Hamiltonian and wave function describing such a 
system are given by 



where 
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¥>(?"i; r 2) is the wave function of two carriers, is the form factor of the electron-phonon 
interaction given by Eq. (jSJ). 

Averaging the Hamiltonian ( |T5l) over Q p h, we obtain the effective Hamiltonian 



H 



B 
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exp(iqr 2 )) + F*(n, r 2 )V* (exp(-iqn) + exp(-igr 2 ))] . (18) 

The displacement amplitude F^(fi,f 2 ) in Eq. (|T8l) . which is a generalization of that 
applied above to the polaron problem, can be written in the form 

F q (n, r 2 ) = f g + Sf g (exp(-igri) + exp(-igr 2 )) (19) 

Substituting Eq. ( |T9l into Eq. ( 1T81) . averaging the Hamiltonian ( {TBI) over <^(ri,r 2 ), and 
minimizing the energy E B = (ip\H B \ip) with respect to f$ and g^, we obtain the following 
functional for the total energy of a 2D bipolaron 



E B = K 1 + K 2 + V 12 - ^ ^< 



v ivg 2 [^-% 2 ] 2 
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where 
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W q = J d 2 r 1 d 2 r 2 [exp(igri) + exp(zgr 2 )] <^ 2 ( r i ; r 2) 

Dg = J d 2 r 1 d 2 r 2 [2 + exp(ig(r x - r 2 )) + exp(-ig(ri - r 2 ))}(p 2 (r 1 , r 2 ) 

The constant of the Coulomb interaction U = e 2 /£oo can be written as U = 2fts/a; 
where /3 — 1/(1 — 77), r\ — Soo/eq. The latter sum in Eq. (1201) . just as in the polaron 
problem, is the correction term arising from the electron-phonon correlation caused by 
the second term in Eq. (fT9l) . The functional (TSUI) without this term (g^ = 0) determines 
the energy of the strong-coupling adiabatic bipolaron. For calculation of the bipolaron 
energy in the strong-coupling limit, we choose the following two trial wave functions: 

¥>(ri,r 2 ) = Nexp(-5 2 (r 2 + r 2 2 ))[l + 7 <* 2 r 2 2 ], 7 > (21) 

<p(r!,r 2 ) = Arexp(-5 2 (r 2 + r 2 ))[l-6exp(- 7 5 2 r 2 2 )], 7 > 0, b < 1(22) 

The bipolaron is stabile when the bipolaron energy is lower than twice the polaron 
energy, i.e., the bipolaron stability regions is determined by the inequality E B < 2E P . 
At a —>■ 00, the bipolaron energy calculated with the use of the trial wave function 
(T2~T]) is approximately 2.5 times larger than the polaron energy E p = —(ir/8)a 2 , 
i.e. E B /2E P m 1.25. While the value of E p calculated in the strong-coupling limit 
(a — > 00) with the use of the trial wave function (122]) is equal to E B ~ 2.58£ , p , so that 
E B /2E P ~ 1.29. In the intermediate coupling regime, the energy of the 2D bipolaron is 
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Figure 2. The 2D bipolaron energy (solid curve), calculated with using the Eq. (120I) . 
and the double 2D polaron energy (dashed curve), calculated with using the Eq. ([9]), 
as a function of a and in units of hio D . 

calculated by using the trial wave function ( 1221) . The results for the limit case (3 = 1 
(f] = 0) are shown in Fig. 2. In this approximation the 2D bipolaron is stable for 
a > a c = 2.94 which is comparable with the 2D results obtained by using the Feynman 
path integral method [9]. With increasing a, the bipolaron stability region is enlarged. 
For example, at a = 5 and 10 the values of the parameter (3 corresponding to the 
bipolaron formation region lie in the intervals 1 < (3 < 1.201 and 1 < (3 < 1.209, 
respectively. For the cases a = 5 and 10 the 2D bipolarons are stable at 77 < r\ c = 0.167 
and rj < r\ c = 0.173, respectively. 

4. Conclusions 

In this work we have studied the formation of intermediate-coupling optical polarons 
and bipolarons in 2D systems. The total energies of the 2D polaron and bipolaron are 
calculated by using the Buimistrov-Pekar method of canonical transformations. It is 
shown that this method allows us to take into account the effects of the electron-phonon 
correlation on the energies of the 2D (bi)polaron for any values of the electron-phonon 
coupling constant a and to derive the new expression for the polaron energy and limits 
for strong and weak coupling. In the strong-coupling limit the Buimistrov-Pekar method, 
the Feynman path integral method and the LLP unitary transformation method give 
the same result for the polaron energy. Our numerical results show that below a critical 
value a c , the charge carriers are quasi-free (i.e. delocalized) electrons or holes, and only 
at a > a c ~ 2.94 (bi)polarons exist in 2D polar materials. The above critical values a c 
and rj c are also close to the ones obtained by using the Feynman path integral method. 
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Note that for calculation of the energies of 2D polaron and bipolaron were used 
simple trial wave functions. Therefore, the values of the (bi)polaron energies, a c and 
r] c obtained by these wave functions are approximate ones. For obtaining more precise 
values of these parameters, it is necessary to use more flexible wave functions, such as, 
for example, a sum of Gussians [2lj . 
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1. Introduction 

As is well known, charge carriers (electrons and holes) is ionic crystals interact with 
polar optical phonons and they are self-trapped at their sufficiently strong coupling to 
the phonons with the formation of optical polarons (called as the Pekar or Frohlich 
polarons) [H EJ [3] . One distinguishes three distinct regimes of electron-phonon coupling 
[2]: (i) the weak-coupling regime describes the correlated motions of the lattice atoms 
and the quasi-free charge carriers which remain in their initial extended state, (ii) the 
intermediate-coupling regime characterizes the self-trapping of a charge carrier which 
is bound within a potential well produced by the polarization of the lattice in the 
presence of the carrier and follows the atomic motions, and (iii) the strong-coupling 
regime describes the other condition of self-trapping under which the lattice atoms 
no longer follow the charge carrier motion and the self-trapping of carriers is usually 
treated within the adiabatic approximation (i.e. lattice atoms remain at their fixed 
positions). Under certain conditions, two charge carriers interacting with the lattice 
vibrations and with each other can form a bound state of two carriers in polar materials 
within a common self-trapping well. Since the attractive interaction of the electron- 
phonon coupling in these systems is strong enough to overcome the Coulomb repulsion 
between two carriers. The self-trapped state of the pair of charge carriers is termed 
a bipolaron. In the last few decades, the bipolaron problem has been a focus of 
attention due to its importance in semiconductor technology and in the bipolaronic 
mechanism of superconductivity (see Refs. [H 13, EJ d El [9], UTIl HH H2])- After the 
discovery of the layered high-T c cuprate superconductors, the study of bipolarons has 
attracted the revived interest because some researchers believe that the bipolaron is 
one of the possible candidate for explaining high-T c superconductivity. The mechanism 
of high-T c superconductivity of large bipolarons was proposed by Emin and Hillary 
[7j. This mechanism is based on the Bose-Einstein condensation of bipolarons as 
discussed by Schafroth [13]. Another mechanism of high-T c superconductivity driven 
by the superfluid single particle and pair condensation of large bipolarons and polaron 
Cooper pairs was proposed in [TJ1E5]. The possibility of such a novel superconductivity 
depends on the existence of polaron and bipolarons in the superconducting materials. 
The above electron-phonon coupling regimes are characterized by the dimensionless 
Frohlich coupling constant a. The possible ranges of the values of a characterizing 
the weak, intermediate and strong coupling regimes depend on the type of self-trapped 
charge carriers and the dimensionality of the system. The polaron and bipolaron ground 
state in the strong coupling limit (i.e. in the adiabatic approximation) have been studied 
by many authors for 3D and 2D systems 0, [7J El EJ EH1 EH ES Ej- However, with 
decreasing a, it is necessary to take into account the electron-phonon correlation which 
reduces the energies of the polaron and bipolaron in comparison with those of the 
adiabatic polaron and bipolaron. Therefore, a quantitative treatment of the (bi)polaron 
problem in the intermediate-coupling regime was also a subject of numerous studies in 
the past three decades (see Refs [61 [8j Ell ESI HI EH]). The ground state energies 
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of the polaron and bipolaron have been calculated by using several approximations, 
such as the Feynman path integral method [HI [9j [101 [HI [20], the Lee- Low-Pines (LLP) 
unitary transformation method [321 [18] and the operator formalism [HI [IE]. According to 
different intermediate coupling treatments, the formation of the polaron and bipolaron 
becomes possible only if the values of a are greater than the certain critical values a c . 
The previous calculations have shown that the (bi)polaron is created more easily in 
2D systems than in 3D ones and the stability region for bipolaron formation is much 
broader in 2D compared with 3D case. So far, calculations of the ground-state 

energies of the (bi)polaron, the values of a c and the stability region of the bipolaron 
are not conclusive. Some of the estimated values of a c differ greatly in magnitude. 
Further, the calculation methods used to study the (bi)polaron problem suffer from the 
certain drawbacks (see Refs. 0,118]). In particular, the Feynman path integral method 
is more accurate for the calculation of the polaron energy pi [161 HB], but it fails to 
describe correctly the bipolaron formation. Moreover, the functional for the bipolaron 
energy has rather a cumbersome form which requires the tedious numerical calculations. 
Therefore, it is expedient to use the another method for the calculation of the polaron 
and bipolaron energies in the intermediate-coupling regime. One of such methods is 
the Buimistrov-Pekar method |21j. In this approach, the electron-phonon correlation is 
taken into account through the displacement amplitude F$(r) of the form: 

F <i(r) = k + 9qexp (-iqr) , ( 1 ) 

where f$ and are the variational parameters, which are determined from the condition 
for minima of the total energy, f is the electron coordinate, if is the phonon wave vector. 

The second term in Eq. ([Q) takes into account the electron-phonon correlation 
effect that leads to the reduction of the energy of a 3D adiabatic polaron when a > a c . 
At a < a c the present method leads to the delocalized state of a polaron, and the 
total energy of the polaron is given by E p = afru (where u Q is the longitudinal optical 
(LO) phonon frequency). In the Feynman path integral variational approach [20], the 
derealization of the polaron at a finite value of a does not occur and the lowest polaron 
energy is obtained. However, at a > a c the Buimistrov-Pekar method taking into 
account an important part of the electron-phonon correlation gives also reasonable 
results for the energy of the polaron. The main advantage with this method is that it is 
simple and does not require tedious numerical calculations. Recently, the Buimistrov- 
Pekar method was applied to study both the free bipolaron and the bound bipolaron in 
3D systems [22l l23l [24] . The existence of polarons and bipolarons in high-T c cuprates 
and other materials has been indicated by several experiments [251 [261 [271 1281 129] . The 
Cii0 2 -based layered high-T c materials are believed to be quasi-two-dimensional systems 
and in an intermediate electron-phonon coupling regime. So far, the problems of the 
2D polaron and bipolaron are not studied sufficiently and the possibility of formation 
of intermediate-coupling polarons and bipolarons in 2D systems within the Buimistrov- 
Pekar formalism is not explored. The aim of the present paper is to study the formation 
of the intermediate coupling optical polaron and bipolaron in 2D systems by using 
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the Buimistrov-Pekar method and to expose the important features of this method. 
We calculate the ground-state energies of the intermediate coupling 2D (bi)polarons, 
the critical values of a, the parameter of the Coulomb repulsion between two carriers 
(3 = 1/(1 — Eoo/eq) (where £ ro and e are the high frequency and static dielectric 
constants) for the formation of 2D polaron and bipolaron. We discuss the obtained 
results and compare them with previous ones obtained by using the Feynman path 
integral method and the modified LLP unitary transformation method [16, 30J. 



2. Formation of intermediate-coupling polarons 



The Hamiltonian and variational wave function describing the interacting system of 
electron (or hole) and LO phonons can be written as [221 El 

h 2 



H„ 



2m 



■ A + E \- V i b i ex P(^ r ) + v q *K exp(-igr)] + ^obp q , (2) 



and 



where 



* = % h ip{r) = U\0)ip{r) 



U = exp 



\V a 



{F q {r)b + q - F q *(r)b q ) 
- g 

2Tthjj l 



U*U = 1, (0|0) = 1, 



-a, 



I, 



h 



a 



2lTUJ l 



1 



1 

CO 



(3) 



(4) 



(5) 



huo L 2 q V 2m*uj 

m* is the effective mass of the carrier before polaron formation, Vg- is the 2D form factor 
of the electron-phonon interaction, L 2 is the size (or surface) of the 2D system, b~t{b q ) 
is the creation (annihilation) operator of the LO phonon with the wave vector q and 
the frequency u , <p{r) is the one-electron wave function represents the phonon part 
of the wave function, |0) is the unperturbed zero-phonon state satisfying b^O) = and 
(0|0) = 1. Averaging the Hamiltonian (JSJ) over $ p ^, we have 



h 2 



2m' 



_5_ 

2m* 



\VF q (r)\ 2 + hw \F g (r)\ 2 + 



V q F q (r) exp(iqr) + V q *F*(r) exp(-iqr)\ (6) 

After substituting Eq.Q into Eq.(j6]), averaging the Himiltonian (jBJ) over (p(r), and 
minimizing the energy E p = (tp\H p \tp) with respect to fg- and g$, we obtain the following 
functional for the total energy of a 2D polaron 



E p = K 



E 



fUjJn 



-w 



E 



\V a \ 



W 2 } 2 



hu 1-W 2 + l 2 Q q 2 



(7) 



where 



K 



h 2 

2m* 



d 2 r 



dip 
dr 
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The latter sum in Eq. j7j) is the correction term arising from the electron-phonon 
correlation caused by the second term in Eq.flT]). The functional (J7J) without this term 
(g$ = 0) determines the energy of the strong-coupling adiabatic polaron. For simplifying 
the calculations, we choose the electron wave function in the form 

<p(r) = Nexp(-5 2 r 2 ), nN 2 = 25 2 J W 2 = exp(-q 2 /A5 2 ) (8) 

Then we obtain the following expression for the energy (in units of hw ) of the 
intermediate coupling 2D polaron: 



o 



r- I [1 - exp(-t 2 )] 2 

E p = 2a* - - 2»a j dt 1 _ eM _ t2)+ ^ H2 , (9) 



where \i = l 8. 

In the weak and strong coupling limits, we obtain from Eq. ([9l) 
~ol for a < a c = 2.94 

E p={ ,2 r ( 10 ) 
1 -far for a — >■ oo 

The Feynman path integral method and the LLP unitary transformation method 
applied to the polaron problem, give the following expressions for the energies of the 
polarons [16]: 



oo 



p 2v 2^2} /D(t) 



o 



w 2 v 2 — w 2 



Dit) = — t + -^s-[l - exp(-trt)] (11) 



and 



oo 

ar A _ f J exp[-(l - j)V/X] 



*r=j-j* ' L fr# - ■ (12) 



o 



where if, w and A, 7 are the respective variational parameters. 

The energies of the 2D polarons in the weak and strong coupling limits follow 
directly from Eqs. (11) and (12): 



■\a — ^a 2 for a ^ 



8 [ cr for a — ► 00 



-a for a < a r = 3.62 



8* 
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Figure 1. The energies of the 2D polarons (in units of huj ) as a function of a for the 
strong coupling (clotted curve) and weak coupling (dashed line) limits, calculated by 
using Eq. (|10[) , for the intermediate coupling limit (solid curve), calculated by using 
Eq.©. Thin solid curve calculated by the Feynman path integral method, Eq. (fTTj) . 



These results can be compared with Eqs. ([9]) and (fit)]) . As can be seen from Eqs. (JTDjl . 
({TBI and (1141) . all three methods in the strong-coupling limit give the same result. For 
compairing, in Fig.l, the energies of 2D polarons (in units of ftw ) obtained from the 
Eqs. ([9]) and ( fTTl) are plotted as a function of a. As seen from Fig. 1, the energy of 
the 2D polaron calculated by taking into account the electron-phonon correlation within 
the Buimistrov-Pekar approximation is lower than that of the 2D adiabatic polaron. At 
a > a c ~ 2.94, a bound state of a carrier exists in 2D systems and the polaron remains 
localized in the 2D polarization well. 



3. Formation of the intermediate-coupling bipolarons 



In the bipolaron problem the total energy of the interacting system of two carriers 
and LO phonons is determined. The Hamiltonian and wave function describing such a 
system are given by 



where 



H 



B 



*-2 2 

(Ai + A 2 ) + + Y^i V M e Miqri) + exp(zgr 2 )) + 



2m? ' £oo?T2 

q 

y 9 *6*(exp(-igr!) + exp(-iqr 2 ))} + ^ ^ob^h 

q 

^ = % h (f{n,r 2 ) = U\0)ip(ri,r 2 ), 



U = exp 



J2(F q (n,r 2 )b+-F;( ri ,r 2 )b g ) 



(15) 
(16) 

(17) 



^(Tj r 2) is the wave function of two carriers, V$ is the form factor of the electron-phonon 
interaction given by Eq. 
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Averaging the Hamiltonian (11 5p over <£> p /i, we obtain the effective Hamiltonian 
H B = -^(A, + A 2 ) + -f- + £ (|V I F t (r.,^)|' + 

|V 2 F g (n,r 2 )| 2 ) +^ |F 9 (r 1 ,r 2 )| 2 + F g (r 1 ,r 2 )\/ g (exp(igr 1 ) + 

exp(igr 2 )) + F*(n, r 2 )V*(exp(—iqr 1 ) + exp(-zgr 2 ))] . (18) 

The displacement amplitude i^(r*i,r 2 ) in Eq. ( fTBI) . which is a generalization of that 
applied above to the polaron problem, can be written in the form 

Fgin, r 2 ) = f q + g q (exp(-iqn) + exp(-zgr 2 )) (19) 

Substituting Eq. ffT9]) into Eq. (TT81) . averaging the Hamiltonian ffT8]) over yj(ri,r 2 ), and 
minimizing the energy Eb = {^\Hb\^) with respect to fq- and gr^, we obtain the following 
functional for the total energy of a 2D bipolaron 

WX„, 2 v-IKI 2 P a -w|' 2 



9 9 



where 



^12 



<9ri i2 

W g = J d 2 rid 2 r 2 [exp(zgri) + exp(zgr 2 )] V9 2 (ri, r 2 ) 

D q = J d 2 r 1 d 2 r 2 [2 + exp(ig(r x - r 2 )) + exp(-zg(ri - r 2 ))]f 2 (r 1 , r 2 ) 

The constant of the Coulomb interaction U = e 2 je^ can be written as U = 2 / 9- v / a; 
where (3 = 1/(1 — 77), 77 = £00/^0- The latter sum in Eq. (1201) . just as in the polaron 
problem, is the correction term arising from the electron-phonon correlation caused by 
the second term in Eq. ( 1T91) . The functional (1201) without this term (g^ = 0) determines 
the energy of the strong-coupling adiabatic bipolaron. For calculation of the bipolaron 
energy in the strong-coupling limit, we choose the following two trial wave functions: 

V (r u r 2 )=Nexp(-5 2 (r 2 1 +r 2 2 ))[l + 1 5 2 r 2 12 ], 7 > (21) 

(p(n,r 2 ) = ATexp(-(5 2 (r 2 + r 2 ))[l -6exp(- 7 (5 2 r 2 2 )], 7 > 0, b < 1(22) 

The bipolaron is stabile when the bipolaron energy is lower than twice the polaron 
energy, i.e., the bipolaron stability regions is determined by the inequality Eb < 2E P . 
At a — > 00, the bipolaron energy calculated with the use of the trial wave function 
(f2~TT) is approximately 2.5 times larger than the polaron energy E p = — (7r/8)a 2 , 
i.e. Eb/2E p w 1.25. While the value of E p calculated in the strong-coupling limit 
(a — > 00) with the use of the trial wave function (122]) is equal to Eb — 2.58E p , so that 
Eb/2E p ~ 1.29. In the intermediate coupling regime, the energy of the 2D bipolaron is 
calculated by using the trial wave function (!22j) . The results for the limit case f3 = 1 
(77 = 0) are shown in Fig. 2. In this approximation the 2D bipolaron is stable for 
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Figure 2. The 2D bipolaron energy (solid curve), calculated with using the Eq. 120 j) . 
and the double 2D polaron energy (dashed curve), calculated with using the Eq. ((9]), 
as a function of a and in units of hu . 

a > a c = 2.94 which is comparable with the 2D results obtained by using the Feynman 
path integral method [9]. With increasing a, the bipolaron stability region is enlarged. 
For example, at a = 5 and 10 the values of the parameter (3 corresponding to the 
bipolaron formation region lie in the intervals 1 < (3 < 1.201 and 1 < (3 < 1.209, 
respectively. For the cases a = 5 and 10 the 2D bipolarons are stable at rj < r\ c = 0.167 
and rj < r\ c = 0.173, respectively. 

4. Conclusions 

In this work we have studied the formation of intermediate-coupling optical polarons 
and bipolarons in 2D systems. The total energies of the 2D polaron and bipolaron are 
calculated by using the Buimistrov-Pekar method of canonical transformations. It is 
shown that this method allows us to take into account the effects of the electron-phonon 
correlation on the energies of the 2D (bi)polaron for any values of the electron-phonon 
coupling constant a and to derive the new expression for the polaron energy and limits 
for strong and weak coupling. In the strong-coupling limit the Buimistrov-Pekar method, 
the Feynman path integral method and the LLP unitary transformation method give 
the same result for the polaron energy. Our numerical results show that below a critical 
value a c , the charge carriers are quasi-free (i.e. delocalized) electrons or holes, and only 
at a > a c ~ 2.94 (bi)polarons exist in 2D polar materials. The above critical values a c 
and r] c are also close to the ones obtained by using the Feynman path integral method. 

Note that for calculation of the energies of 2D polaron and bipolaron were used 
simple trial wave functions. Therefore, the values of the (bi) polaron energies, a c and 
rj c obtained by these wave functions are approximate ones. For obtaining more precise 
values of these parameters, it is necessary to use more flexible wave functions, such as, 
for example, a sum of Gussians [24] . 
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